Energy Dissipation in Driven Granular Matter in the Absence of Gravity 
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We experimentally investigate the energy dissipation rate in sinusoidally driven boxes which are partly filled 
by granular material under conditions of weightlessness. We identify two different modes of granular dynamics, 
depending on the amplitude of driving, A. For intense forcing, A > Aq, the material is found in the collect- 
and-collide regime where the center of mass of the granulate moves synchronously with the driven container 
while for weak forcing, A < Aq, the granular material exhibits gas-like behavior. Both regimes correspond 
to different dissipation mechanisms, leading to different scaling with amplitude and frequency of the excitation 
and with the mass of the granulate. For the collect-and-collide regime, we explain the dependence on frequency 
and amplitude of the excitation by means of an effective one-particle model. For both regimes, the results may 
be collapsed to a single curve characterizing the physics of granular dampers. 

PACS numbers: 45.70.-n, 46.40.-f, *43.20.Tb, 07.10.-h, 46.40.-f 



Introduction. When containers filled with granular material 
are subjected to vibrational motion, a number of interesting 
phenomena is observed. The most prominent among them are 
self-organized convection flows and various segregation ef- 
fects which are reported in a large body of literature, e.g. 03- 
3|. Most, if not all, of these effects found in driven granular 
systems are influenced by gravity, therefore, to study these ef- 
fects in absence of gravity, experiments have been performed 
under conditions of weightlessness in parabolic flights, drop 
towers, and sounding rockets. Examples for such investiga- 
tions concern shear flow [4|, cooling and clustering in dilute 
systems ||5}{8]], violations of the energy equipartition j9 ffl"Tl . 
the propagation of sound [ 12 1 and segregation fl3l . 

In this letter, we experimentally address the mechanisms of 
energy dissipation in vibrationally driven granular systems in 
the absence of gravity. 

Dissipation mechanisms belong to the fundamental proper- 
ties of matter and are of interest per se. But besides scientific 
curiosity, the investigation of dissipation in granular systems 
is of technological interest to dampen unwanted vibration us- 
ing devices which are referred to as granular dampers which 
are containers or cavities filled with granular material. When 
the container is subjected to oscillatory motion, the grains in- 
side the cavity collide inelastically with one another and with 
the confinement and, thus, dissipate mechanical energy of the 
vibration. In contrast to conventional dashpot dampers, granu- 
lar dampers do not rely on a fixed anchor as an impulse reser- 
voir. They show only a minute dependence on temperature, 
and since such dampers can be sealed off hermetically, they 
are predestined for long term use in harsh environment with 
extreme temperatures and/or high pressures. For example, us- 
ing granular dampers, turbine blades can be kept from oscil- 
lating and in medical tools where sterilization is mandatory 
they can dampen the vibration of handles lfT4HT71l . 

Extensive research has been performed on understanding 
the behavior of granular dampers in specific applications and 
to determine the influence of material and construction param- 
eters like the size of the container or cavity, the number of par- 
ticles, the mass of the filling, and the clearance, e.g. Ifl8lj26l . 



Recent simulations show the flow of the granular material in- 
side such a cavity under normal gravity for a vertically shaken 
system Il24ll25l 

From all of these investigations it may be concluded that 
the dissipation properties of a granular damper depend in a 
non-trivial way on many variables like the amplitude of the 
oscillation, the clearance, the frequency as well as the mass 
and type of particles l27ll28l . By now, however, these depen- 
dencies are poorly understood. 

Besides the lack of understanding the physics of such sys- 
tems, from a technological point of view there is a need for a 
master design curve, in order to be able to predict the damping 
properties of a granular damper in a certain application. 

As shown in the comprehensive work by Yang 12T1 for ver- 
tical excitation in gravity, only for strong forcing (high exci- 
tation energy), it is possible to collapse the dissipated power 
as a function of the amplitude for various frequencies and, 
thus, to deduce an empirical master design curve. For less in- 
tense driving it was shown by means of numerical simulations 
||29ll that to a large extend the dissipative properties of granular 
dampers are influenced by gravity, in particular when gravity 
exceeds the acceleration due to the external forcing. 

To isolate the response of driven granular systems to ex- 
ternal excitation from the disturbing effect of gravity, in this 
letter for the first time we address the energy dissipation of 
granular dampers under conditions of weightlessness during 
parabolic flights. We explain the dependence on frequency 
and amplitude of the excitation by means of an effective one- 
particle model. Without using any adjustable parameter the 
results may be collapsed to a single curve characterizing the 
physics of granular dampers. 

Experiment. The experimental setup is sketched in Fig. [T] 
The sample box is mounted on a strain gauge which in turn 
is attached to a carrier moving on a linear bearing. A gear 
belt connects the carrier with a computer controlled stepper 
motor which drives the carrier to perform sinusoidal oscilla- 
tions of adjustable angular frequency, uj, and amplitude, A. 
The time-dependent position of the container was measured 
using Hall-effect based position encoders with a resolution of 
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20 ^m and 10 kHz sample rate to check that the deviation of 
the container's motion from the set sinusoidal oscillation is 
negligibly small, typically the spurious free dynamic range is 
56 dB. The experiment was monitored by a high-speed camera 
at a frame rate of 240 fps. The entire setup was built up twice 
such that the containers move in opposite direction to cancel 
the vibrations transmitted to the external mounting structure. 



motor 




box 
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linear bearing Q ^ 



FIG. 1 : Sketch of the experiment. For explanation see the text. 

The strain gauge delivered signals proportional to the forces 
in the direction of driving, and the two directions perpendicu- 
lar to it. Only the first one is relevant for our experiment; we 
checked that the forces in the other directions are negligibly 
small as compared to it, i.e., the side walls confine the granular 
material while the main transfer of momentum occurs parallel 
to the direction of motion. 

The samples consist of polycarbonate boxes (wall thickness 
4 mm) partially filled by different amounts of steel beads (di- 
ameter 4 mm, material density 7.8 g/cm 3 , Young's modulus 
203.5 GPa). The number and the total mass of particles within 
each sample is given by N and m, respectively. The clearance, 
L g , is the difference between box length and the thickness of 
the packed layer of particles in the box. It can be obtained by 
computing the volume occupied by particles in random close 
packing lf30ll at volume fraction 64 %. Table[I]summarizes the 
characteristics of our samples. 

TABLE I: Table of Samples. 



Sample No. 


BoxLxWxH(mm 3 ) 


m (g) 


TV 


L g (mm) 


1 


100 x 50 x 50 


126.3 


473 


89.4 


2 


50 x 50 x 50 


135.3 


507 


38.7 


3 


50 x 50 x 50 


71.0 


266 


44.1 


4 


100 x 50 x 50 


63.8 


239 


94.7 



era. For the results reported here we use only the data obtained 
in the stationary state. 

Regimes of Dynamical Behavior. Analyzing the high-speed 
video recordings we can identify two different regimes of dy- 
namical behavior, see Fig. [2] For large amplitudes of the vi- 
bration, the damper operates in the collect-and-collide regime, 
that is, during the inward stroke all the material is "collected" 
and accumulates as a relatively densely packed layer at the 
wall of the container. After passing the phase of maximal ve- 
locity, the box decelerates and the layer of particles leaves the 
wall collectively. When the bulk of particles impacts the op- 
posite wall of the container, a large part of the kinetic energy 
is dissipated by inelastic collisions. The amount of energy 
dissipated depends on the relative velocity between the parti- 
cles and the wall at the time of the impact, determined by the 
amplitude and frequency of the vibration and the filling ratio 
of the container. The collect-and-collide regime was theoreti- 
cally and experimentally investigated in [31 1 and identified as 
the regime of most efficient damping. It was confirmed also by 
numerical MD-simulations ||32l [33l and identified as one out 
of four different regimes of dynamical behavior of vibrated 
granulate in microgravity. 




To exclude the influence of gravity, the experiment was per- 
formed during a parabolic flight allowing for stable micro- 
gravity condition (0 ± 0.05) g for time intervals of about 22 
seconds which determines the duration of each single mea- 
surement, where amplitude and frequency of the excitation 
were fixed. The data resulting from the strain gauges and 
the position sensors were sampled simultaneously at a rate of 
10 kHz and stored for later evaluation of the dissipated energy. 
About four seconds after the onset of microgravity the experi- 
ment had entered the stationary state which could be identified 
by both the rate of dissipation deduced from the measurement 
of the driving force and the recordings of the high-speed cam- 



FIG. 2: Snapshots from the high-speed video recordings of sample 4 
illustrating the two distinct regimes of excitation: collect-and-collide 
regime at A — 50 mm (left column) and a gas like state at A = 
2.5 mm (right column). Each column shows the box at a phase range 
from (top) to n (bottom). 

For small amplitudes we observe a gaseous state where 
only a small fraction of the particles interact with the oscillat- 
ing walls during one oscillation period. In the gaseous state, 
the collisions of the particles with the driving walls are just 
sufficient to balance the energy loss according to dissipative 
particle-particle collisions in the bulk of the material. Here, 
the dissipation rate is smaller than in the collect-and-collide 
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regime fl3~Tl . 

Energy Dissipation Rate. To obtain the energy dissipated 
by the granulate during one period, T = 2ir/us, of the sinu- 
soidal driving, x = A sin(wi), we integrate the product of the 
measured force, F(t), and velocity, x(t), over one period of 
oscillation: 



E dlss = J x(t)F(t) dt. 



(1) 



The maximum energy that can be dissipated during one cy- 
cle in the system is given by: 



12, ,2 



(2) 



This is the case if all particles collide inelastically with the 
wall at maximum relative velocity. In the following, E max is 
used for normalization. 

We measured the dissipated energy per period for the fol- 
lowing ranges of frequency and amplitude: Sample 1 and 4 
were shaken at 1 Hz, 2 Hz, 4 Hz while samples 2 and 3 were 
shaken from 1 Hz to 5 Hz in 1 Hz increments. For each setup, 
Fig. [5] shows -Bdiss/^max versus the amplitude of the oscilla- 
tion. 

Let us first consider the gas-like state observed for small 
amplitude, A < A . In this regime, we expect the dissipated 
energy to be proportional to the number of particles collid- 
ing with the wall. If we assume a monodisperse system with 
homogeneous density this number is determined by the the 
volume swept by the container's wall. We further assume the 
characteristic velocity of the particles to scale with the veloc- 
ity of driving, Alu, and the particles hitting the wall at random 
phases, due to their disordered motion and arrive at: 
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FIG. 3: (color online) Normalized dissipated energy per period of 
external vibration. Symbols: Experimental data. O : 1 H z > □ : 
2 Hz; <) : 3 Hz, < : 4 Hz; t> : 5 Hz. Lines: Solution of the impact 
model, valid for A > Aq (solid - numerical, dashed - analytical, Eq. 
l(9](). Dotted lines: Dissipation rate for the gaseous regime (A < Ao, 
see Eq. {3)). Inset: same data (only gas regime) but normalized to 
E mm /m (see Eq. Q). The error bars for the gas regime are shown in 
the inset. For all other measurements, the errors are about the size of 
the symbols. The threshold amplitude, Aq, (vertical lines) obtained 
from the model (see Eq. agrees with the experimental data. 
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Note that particle-particle collisions in the bulk of the ma- 
terial contribute only indirectly to £2 since such collisions 
do not transfer momentum to the container. 

Equation |3]l was developed under the assumption of a ho- 
mogenous density distribution. This however may not always 
hold true. Unlike molecular gases, heated granular gases are 
not homogeneous but density increases in a non-linear way 
with distance from the driving wall ll34l to form regions of 
enhanced density (clusters) far away from the wall. Follow- 
ing the arguments of ||35l the number of particle wall collision 
depends only weakly on the total mass of particles in the sys- 
tem. Consequently, for the limit of no dependence on the total 
mass we may write 
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shown in the inset of Fig. [3^. 

For the cases described by Eqs. ( 3]4 1, the experimental data 
shown in Fig. [3] collapse despite the fact that the data points 



shown for a certain amplitude correspond to different frequen- 
cies of driving. That is, the dissipation rate is independent of 
frequency in agreement with the scalings, Eq. (|4]) for Fig. [3]a) 
and Eq. ^ for Fig. |3]b). 

For large amplitudes, A > Ao, when the granulate is in 
the collect-and-collide regime, the experimental data collapse 
when scaled with E max . Moreover, also here the scaled data 
belonging to the same amplitude are independent of the fre- 
quency. To quantitatively explain the data shown in Fig. [3] and 
to understand its independence of frequency for A > Ao, we 
take a closer look at the dynamics of the collect-and-collide 
regime, where nearly all particles are collected by the inward 
stroke and leave the wall collectively at t = with veloc- 
ity %micie = Aw. We define t c as the time when the parti- 
cles collide with the opposing wall where they will adopt its 
instantaneous velocity w wa u = Aiu cos(wi c ). Note that this 
corresponds to the motion of a quasi-particle interacting per- 
fectly inelastically with the container walls. For a justification 
of this model see ll24l I3T1 l36l . The amount of kinetic energy 
lost per period depends on the difference of the velocity of the 
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quasi-particle and the wall: 

^dfss = m ( w paiticle - Wwall) 2 ■ 

Expressed in terms of t c and E max we obtain: 
^ ss = 7[l-cos(a;t c )] 2 S max 



(5) 



(6) 



The time t c is implicitly given by the distance the bulk of 
particles has traveled and the harmonic motion of the box, 



^waii^c = Aujt c = Asm(cjt c ) + L g . 



(7) 



Equation Q can be solved numerically for uit c to obtain the 
dissipated energy per period, E%? ss /E max , via Eq. (|6]l (see solid 
lines in Fig. [3J. 

Alternatively, we can obtain an approximate value for uit c 
by a first order expansion of Eq. |7]i around u)t c = 7r: 



Ujt r 
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Inserting this solution into Eq. (|6]l we find 
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Figure [5] compares the relative dissipated energy per oscilla- 
tion period as obtained in experiments (symbols) with the nu- 



merical solution of the collect-and-collide model, Eqs. (5]7l, 
(solid lines) and its approximate analytical solution, Eq. (|9]l 
(dashed lines), resulting in good agreement. Also in agree- 
ment with the experimental data, both Eq. |9) and Eq. (|6]l 
with u)t c from Eq. |7]i and Eqs. ( |3|4| i are independent of the 
frequency, uj, which explains the collapse of the data for dif- 
ferent frequencies. Note that the model, Eqs. ( |6|7| >, and the 
approximate solution, Eq. (|9), do not contain any adjustable 
parameters. 

Furthermore, this model provides an explanation of the 
threshold A separating the gaseous state from the collect- 
and-collide regime: The bulk of particles leaves the wall on 
its inward stroke at time t = (instant of maximal velocity). 
If it arrives at the opposite wall at a time where the wall is 
accelerating away from it, that is it < u)t c < 2ir, it will not 
get collected by the wall but mainly scattered. This scatter- 
ing inevitably leads to desynchronization of collective particle 
motion and the collect-and-collide mode breaks down. The 
threshold can be obtained from Eq. ([H} with ujt c = n: 



A 



(10) 



Interestingly, at the edge of stability of the collect-and-collide 
regime, A = Aq, that is, uit c — ir, from Eq. |6) we obtain 
maximal efficiency in damping. The threshold amplitude, Ao, 
obtained from the model agrees well with the experimental 
data, see Fig. [3](vertical lines). 

Conclusion. On Earth, the influence of gravity on granu- 
lar dynamics can only be neglected for intense driving, where 



A 2 lu 2 g. To study the response of granular dampers to ex- 
ternal excitation isolated from the disturbing effect of gravity, 
we investigated the energy dissipation of granular matter in 
the absence of gravity when subjected to sinusoidal motion, 
x(t) = A sin(cji). 

Depending on the amplitude of the vibration, we observe 
two qualitatively different modes: For small amplitude, A < 
Aq « Lg/ir, the granular material behaves gas-like while for 
larger amplitudes, A > Aq, we observe a collect-and-collide 
regime where the center of mass of the granulate moves syn- 
chronously with the driven container. In this regime, we de- 
scribe the system by an effective one-particle model, Eq. |6|, 
which can be solved numerically or using an analytical ap- 
proximation. Despite of the model's simplicity, we find re- 
markably good agreement with the experimental data. 

In the gas regime, the energy dissipation rate increases lin- 
early with the amplitude A. For large clearance, L g , the en- 
ergy dissipation rate is almost independent of the mass of the 
granulate in agreement with the known density instability in 
granular gases, Eq. Q. Here an addition of particles hardly 
increases the density in the vicinity of the heated wall l35l . 
For smaller clearance, the high-density region is no longer far 
away from the driving wall such that the energy dissipation 
rate increases linearly with the mass, Eq. Q 

In both regimes, gas-like and collect-and-collide, the en- 
ergy dissipation rate was found independent of the frequency, 
in agreement with our model description, Eqs. (4|6l. 

With appropriate scaling, for a fixed container geometry all 
measurements coincide along one universal "master curve", 
independent of filling fraction, mass of the granulate, particle 
properties, particle number, and the parameters of vibration. 
Based on the theory of granular gases for the gas state and 
on the effective one-particle model for the collect-and-collide 
state, we derived a model description without any fitting pa- 
rameters which describes the experiment up to a good agree- 
ment. Although not measured, from the model it may be ex- 
pected that the dissipative properties are invariant of the type 
of granulate, the particle diameter and the exact cross section 
of the cavity. 
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